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ABSTRACT 

We show that the photonic analogue of the Korringa-Kohn-Rostocker method is a viable alternative to 
the plane-wave method to analyze the spectrum of electromagnetic waves in a three-dimensional periodic 
dielectric lattice. Firstly, in the case of an fee lattice of homogeneous dielectric spheres, we reproduce the 
main features of the spectrum obtained by the plane wave method, namely that for a sufficiently high 
dielectric contrast a full gap opens in the spectrum between the eights and ninth bands if the dielectric 
constant s s of spheres is lower than the dielectric constant Sb of the background medium. If e s > Eh, 
no gap is found in the spectrum. The maximal value of the relative band-gap width approaches 14% in 
the close-packed case and decreases monotonically as the filling fraction decreases. The lowest dielectric 
contrast Sb/s s for which a full gap opens in the spectrum is determined to be 8.13. Eventually, in the 
case of an fee lattice of coated spheres, we demonstrate that a suitable coating can enhance gap widths 
by as much as 50%. 
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1 Introduction 



Under certain conditions, a gap can open in the spectrum of electromagnetic waves in a dielectric medium, 
independent of the direction of their propagation |^| . Dielectric structures possessing such a photonic 
band gap are promising candidates for various technological applications Q. Moreover, such structures 
offer a new laboratory to study various atomic processes. Indeed, if a gap opens in the spectrum of 
electromagnetic waves, all parameters and characteristics of atom placed in such a medium, such as, for 
example, atomic radius and its spontaneous emission rates are expected to change. 

In order to open such a gap, one considers Maxwell's equations in a dielectric with a spatially periodic 
dielectric function, in full analogy to the Schrodinger equation with a periodic potential [0, ||. In the 
latter case, the spectrum can be classified according to the Bloch momentum k. Energy (frequency) 
levels v n are continuous functions of the Bloch momentum k in the (first) Brillouin zone. We say that 
there is a full gap, or simply a gap between the nth and (n + l)th levels when v n+ i(k) > i/ n {k!) for all 
k and k'. We say that there is a direct gap between the nth and (n + l)th levels when f ra+1 (k) > ^„(k) 
for all k. For the Schrodinger equation in one space dimension, the number of gaps is in general infinite 
and the only periodic potential which does not open any gap in the spectrum is a constant potential || . 
However, the situation changes dramatically in two and higher dimensions. One can prove rigorously that 
the number of gaps in the spectrum can only be finite and, if the potential is not strong enough, no gap 
opens in the spectrum Q . If electromagnetic waves are considered, opening a gap in the spectrum is even 
more difficult and it took several years of intensive search to achieve it experimentally for microwaves [^) . 
Note that Maxwell's equations enjoy scale invariance so that, in principle, by scaling all sizes of a given 
structure one can shift a gap theoretically to whatever frequency range. 

So far, the plane-wave method §, @, |, |, has been the main tool to calculate the spectrum 
of electromagnetic waves in three-dimensional dielectric lattices. However, the plane-wave method is 
numerically rather unstable for a setup considered in experiment, namely when the dielectric function is 
piece- wise constant and changes discontinuously j^, O, Oj. The main culprit for this behaviour is the 
Gibbs instability - the dielectric constant is poorly estimated near its spatial discontinuities by a truncated 
Fourier series which wildly oscillates even if more than one thousand of plane waves is retained (see, for 
example, figure 2 of Q). Also, the plane wave expansion becomes impractical if the dielectric constant 
exhibits a large imaginary part. Another approach to calculate the spectrum of electromagnetic waves in 
three-dimensional dielectric lattices uses a discretization of Maxwell's equation inside the primitive cell 
of the lattice pi. However, both methods are difficult to apply in the presence of impurities and to the 



calculation of Green's function. 

In order to have an universal method which can deal with problems of the behaviour of electromagnetic 
waves in a periodic dielectric medium in their full complexity, we have developed and employed a photonic 
analogue of the first principle on-shell multiple-scattering theory (MST) and of the Korringa-Kohn- 
Rostocker (KKR) method |Tq| . The unique feature of the on-shell MST is that, for nonoverlapping 
(muffin-tin) scatterers Jl7| (the present situation), it disentangles single-scattering and multiple-scattering 
effects (sec [jl8| for a recent discussion). The KKR method [|l6| uses explicitly scattering matrices and 
Green's function which are expanded in the basis of spherical harmonics and the spectrum is determined 
by zeroes of a determinant. For electrons on a Bravais lattice, inclusion of spherical waves with angular 
momentum up to l m ax = 2 already gives result within a few per cent of the exact calculation [ fl6| |. 
Expansion in the basis of spherical harmonics does not mean that scatterers have to be spherically 
symmetric. Indeed, scatterers of arbitrary shape are allowed in which case scattering matrices are simply 
nondiagonal in the angular momentum indices [^0| . The main advantage of the KKR method is that it 
gives directly the total Green's function from which the density of states (DOS) and the so-called local 
density of states can be easily extracted. The local DOS, which is proportional to the imaginary part 
of the total Green's function at the coinciding points in the coordinate space, is an important quantity 
which determines decay of excited states of atoms and molecules embedded in the lattice [jljj. Also, the 
frequency dependence of the dielectric constant can be easily implemented in the formalism. 

The outline of our paper is as follows. In the following section, we show that the photonic KKR 
method is a viable alternative to calculate the photonic band structure by reproducing the main features 
of the spectrum obtained by the plane- wave method. The lowest dielectric contrast £b/e s for which a 
full gap opens in the spectrum is found to be 8.13 and is slightly lower than 8.4 obtained by the plane- 
wave method JhJ. In Section 3 we discuss the case of coated spheres, i.e. spheres made out of several 
spherical shells with different dielectric constant. We demonstrate that already a suitable single coating 
can enhance some of the gap widths by as much as 50%. Our conclusions are summarized in Section 4. 

2 Face-centered cubic lattice of dielectric spheres 

In this section, we shall present the results of our numerical calculation for a face-centered cubic (fee) 
lattice of dielectric spheres with a single sphere per lattice primitive cell. This case is very interesting from 
the experimental point of view, since such dielectric lattices form when silicon matrices, synthetic opals, 
and collodial crystals are used O, [ill IHl. Some of the structures were shown to exhibit the so-called 



stop gap (gap in the spectrum at a fixed direction of the incident light) at optical frequencies [ [ll[ and 
are the natural candidates to achieve a full photonic band gap []12| . 

At the same time, the case of fee lattice of dielectric spheres has been controversial since the first 
experimental results were published [ plf . Results for a sample consisting from polycrystalline AI2O3 
spheres, 6mm in diameter with a microwave refractive index of 3.06 in thermal-compression-molded 
dielectric foam of refractive index 1.01 indicated the presence of a "photonic band gap" in the microwave 
regime pl|| . However, subsequent numerical calculations using the plane wave method |6|, [?J claimed that 
no gap opens in the spectrum and only a pseudo-gap (a sharp drop in the DOS) exists ||. Nevertheless, 
two years later using the plane-wave method, Soziier, Haus, and Inguva [|| did find a full gap for the fee 
lattice of dielectric spheres between the eights and ninth bands. The discrepancy between the results of 
H f?J and H follows from the fact that unlike to the case of electrons, a gap for electromagnetic waves 
opens in an intermediate region and the authors of stopped their calculation just beneath that 

region (see figure 1 in || |7|). Later, the results of Soziier, Haus, and Inguva Q were confirmed by two 
other groups using the plane- wave method §,0. 

The latter deserves some discussion. In the case of electrons, the formation of bands results from 
the broadening of individual atomic levels when the atoms start to feel the presence of each other. The 
largest gap between atomic levels is between the lowest-lying energy levels. Therefore, for a lattice of 
atoms, one expects to find a gap essentially between the first and the second energy band with the gap 
between higher bands scaling down to zero |3j. However, for electromagnetic waves a gap does not open 
between the lowest lying bands but in an intermediate region. This phenomenon can be rather easily 
understood, since for a dielectric scatterer and Maxwell's equations bound states are absent. They are 
replaced by resonances and the above argument for locating the position of a band gap no longer holds. 
Moreover, if the wavelength is small compared to the size of the spheres, one can use geometric optics 
while in the opposite limit of long wavelengths, the Rayleigh approximation applies. In neither case a gap 
opens in the spectrum. Therefore, if a gap is present in the spectrum, it should be in the intermediate 
region between the two limiting cases (see, however, the case of a diamond lattice (Q, figure 2), which is 
a complex lattice). The very same is also expected to apply for the localization of light [22] ]. 

2.1 Results 

Soziier, Haus, and Inguva [^) were well aware of the convergence problems of the plane-wave method and 
they called for the recalculation and confirmation of their results by a more precise method. The latter 
constitutes the first part of our results. Using the photonic KKR method, we were able to confirm the 



plane- wave method result |8j that, in the case of air spheres and for a sufficiently high dielectric contrast 



• a full gap opens between the eights and ninth bands 

• a direct gap opens between the fifth and sixth bands. 

If the dielectric constant of spheres is larger than that of a background medium, no gap opens in the 
spectrum. This situation is realized, for example, if dielectric spheres in air are considered. We did not 
find any compelling explanation for this behaviour. In general, the higher frequency the higher l max is 
to be taken. Taking l max — 1 is sufficient to account for the linear part of the band structure around 
the r point. The intermediate region requires then l max = 3 — 5 and l m ax = 6 is needed to ensure good 
convergence in the range considered. 

In figures land | we present our results for a three-dimensional close-packed fee lattice of air spheres 
in a dielectric medium with the dielectric constant £{, = 9. We choose this configuration for two reasons. 
First, it is sufficiently representative to show the presence of a full gap in the spectrum, and secondly, the 
value of the background dielectric constant Eb = 9 is close to that of rutile (TiC>2) at optical frequencies 
which is used in experiments, figure ^ shows the band structure. Frequency v is plotted in scale-invariant 
units c/irA, where A is the lattice constant^] and c is the speed of light in the vacuum. 

Only a single gap with a middle of gap frequency v — 2.796 and the width Ais = 0.044 opens in the 
spectrum in the range considered. The error is determined from the convergence properties of the KKR 
method. In the close-packed case, the lower gap boundary takes on its maximal value at the W point of the 
Brillouin zone while the upper gap boundary takes on its minimal value at the X point, in agreement with 
the plane- wave calculations (see |24| for the classification of special points of three-dimensional lattices). 
In general, the photonic bands show much more branching than the electronic bands and the actual 
classification of different bands can be quite involved. Group-theoretical classification of eigenmodes in 
three-dimensional photonic lattices is discussed in |25| . 

The presence of the gap in the Eb — 9, e s — 1 case is also transparent from the calculation of the 
DOS per primitive unit cell. The latter was calculated using the Monkhorst-Park integration scheme |p6| . 
Integration over the Brillouin zone started from a mesh of 12 x 12 x 12 uniformly spaced points, which 
was subsequently reduced to 182 points with calculated weigths using the symmetries of the lattice. The 
resulting DOS per primitive cell is plotted in figure ||[ 



1 Noto that A is the side of the conventional unit cell of the cubic lattice, which has four times the volume of a primitive 
fee unit cell, and not the lattice spacing [ pj| . 



Figure ^ shows the band structure for a close-packed foe lattice of air spheres in a dielectric medium 
with the dielectric constant Sb — 2.1609 (rib = 1-47). The latter case corresponds to the experimental 
setup of |llj and is also close to that of fl^ ]. Our calculation shows no gap in the spectrum. Only 
stop gaps are present. In agreement with the experimental observation, the most pronounced stop gap 
is seen between first bands at the L point of the Brillouin zone. For comparison with experiment, in 
Tab. I we give the width Az^ of the stop gap at the L point and the effective refractive index n e // for 
a close-packed fee lattice of air spheres in background media with rib — 1.33, 1.37, 1.47, and 1.6 used in 
recent experiments [ O , |l2]j , together with the case rib = 3 for which the band structure was calculated 
(see figure [l]) . 

TABLE I. The width Az/j, of the stop gap at the L point, effective refractive indices n e // and «-e/f G j and 
strength parameters e r and \& for a close-packed fee lattice of air spheres in different background media. 





rib = 1-33 


rib = 1-37 


n b = 1-47 


rib = 1.6 


rib — 3 




0.145 


0.159 


0.195 


0.236 


0.368 


n eff 


1.084 


1.094 


1.120 


1.153 


1.567 


„MG 
U eff 


1.085 


1.096 


1.122 


1.158 


1.607 


e r 


0.281 


0.313 


0.391 


0.487 


1.140 




-0.376 


-0.410 


-0.485 


-0.566 


-0.935 



The effective refractive index n e ff is determined as the inverse of the slope of the band structure around 
the r point, 

1 di,, 

(1) 



_i 1 du 
n „ — hm , 

e ff k^o c dk ' 



where ui is the angular frequency. In the third row of Tab. I, we show the refractive index nMP calculated 



by the Maxwell-Garnett formula |27 



n eff 



f 2e b + e s + 2f{e s -e b ) 
V 2e b + e s - f(s s - s b ) 



1/2 



(2) 
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where / is the filling fraction (/ = 0.7405 for a close packed fee lattice). In accordance with the plane- wave 
results p8| (see figure 2 there) , n^j^ gives the upper bound on n e ff ■ For < s s the situation is reversed 
and is expected to give the lower bound on n e ff |2S[| . For completeness, we also show parameters 



fel + {l-f)el 



[fe s + (1 - f)e b f 



(3) 



introduced in |b||, and 



introduced in p2| , which should characterize the scattering strength of a dielectric lattice. 

From the experimental point of view, it is interesting to know what is the threshold dielectric contrast 
Smax/^min, where £mox (smin) is bigger (smaller) of the e s and Eb, for which a full gap opens in the spectrum. 
Obviously, this threshold value changes with the radius of spheres and also depends on whether the 
dielectric constant of spheres is larger or smaller than that of the background medium. The precise value 
of the threshold dielectric contrast has been out of reach of the plane- wave calculations Using the 
photonic KKR method, we scanned different configurations between the X and W points of the Brillouin 
zone. For close-packed air spheres, the lower and upper bounds of the full gap are set at the W and X 
points, respectively. For smaller filling fractions, already at / = 0.70, and close to the threshold dielectric 
contrast, the gap width is completely determined by the band structure at the W point. We determined 
the lowest threshold dielectric contrast £b/s s for an fee lattice of dielectric spheres to be 8.13 (e r = 1.096 
and \& = —0.918 in this case). This can for example be realized for the case of close-packed air spheres 
in a background dielectric medium with the dielectric constant £5 = 8.13. In all other cases, i.e., if the 
radius of spheres is lowered, the threshold dielectric contrast is higher. The threshold dielectric constrast 
obtained by the photonic KKR method implies the threshold refractive index constrast 2.8506 which is 
significantly higher than the early theoretical estimate 1.21 by Yablonovitch [jjj and 1.46 by John On 
the other hand, the threshold dielectric constrast is slightly lower than 2.9 obtained by the plane wave 
method @. 

The plot of the relative gap width, which is the gap width divided by the midgap frequency, as a func- 
tion of the the refractive index contrast for different filling fractions is presented in figure |J. The maximal 
value of the relative gap width approaches 14% in the close-packed case and decreases monotonically as 
the filling fraction decreases. The relative gap width as a function of the refractive index contrast shows 
a rapid saturation. For example, in the close-packed case the relative gap width at 5.48 is already 80% 
of its maximal value. 
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3 Face-centered cubic lattice of coated spheres 



After showing in the previous section that the photonic analogue of the KKR method is a viable alternative 
to the plane-wave method to calculate a photonic band structure, we shall proceed by investigating the 
case of coated spheres. A coated sphere is a sphere with the dielectric constant e\ and the radius T\ 
embedded in a larger sphere with the dielectric constant £2 and the radius T2 > r\ , which in turn can be 
embedded in a larger sphere, etc. Let £& be, as above, the background dielectric constant. Our interest 
in this system comes from the fact that there is an intensive experimental activity to produce lattices of 
coated spheres by using them as basic particles in collodial crystals. As in the previous section, due to its 
experimental relevance, we shall investigate only a simple fee lattice (one scatterer per unit cell). It turns 
out that a coating alone does not help to create full gap extending over all the Brillouin zone as long 
as an fee lattice remains a simple lattice. Nevertheless, following a suggestion by A. van Blaaderen, we 
shall show that a suitable coating of spheres forming a dielectric lattice can significantly enhance (by as 
much as 50%) some of the stop gaps. We looked mainly at the so called lowest L-gap width (see |24j for 
the classification of special points of three-dimensional lattices), which corresponds to the (111) crystal 
direction (see for a related theoretical discussion of the L-gap). The reason is that experimental 
techniques make it possible to allow one to grow collodial crystals such that the L direction corresponds 
to normal incidence on the crystal surface. Another reason is that recently pQ |, in the case of air spheres, 
a simple formula has been found for the L-gap width of an fee lattice. For the case of dense spheres 
such an understanding is missing. There is a hope that fee lattices of coated spheres could provide some 
insight into the L-gap behaviour. 

At first it seems that a coating will not have any significant effect on gap widths. One can arrive at 
this conclusion by looking at the effective dielectric constant e c e ^ for a lattice of coated spheres. As we 
have seen in the previous section, the effective dielectric constant is a measure of the scattering strength 
of a dielectric lattice almost up to the first stop gap (Bragg reflection peak). Let us begin with the case 
of a coated sphere with a single coating, in which case we shall denote the dielectric constant and radius 
of the outer shell by e and r s , respectively. Defining x — rf/rf, 

£1 + 2e £ + lEb 

the polarization factor a c of a coated sphere with a single coating is 

a Q + xai(e b + 2e)/(s + 2e b ) 



1 + 2xotiOtQ 
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(6) 



Here < x < 1 and < Oij < 1. The Maxwell-Garnett formula [^7| then gives 

£ J f «£i(l + 2W/(l-/4 (7) 

One can verify that for x = 1, i.e. if the radius of the interior sphere coincides with that of the entire 
sphere, 

El + 2£fa 

i.e. one reproduces the polarization factor of a homogeneous sphere with the dielectric constant e\ in 
the host with the dielectric constant e&. The same result one also recovers in the special case e\ = e. In 
the limit x — ► 0, i.e. if the radius of the interior sphere shrinks to zero, a c — * eta, the polarization factor 
of a homogeneous sphere with the dielectric constant e in the host with the dielectric constant It is 
important to realize that as a function of the parameter x, a c (x) is a monotonic function continuously 
interpolating between c*o and a,\. The latter follows from the fact that the first derivative of da c (x) / dx 
has a constant sign determined by the sign of e — Sb- Since 



d£ eff _ 3M 



> 0, (9) 



da c (1 - fa c ) 2 

also the effective dielectric constant e c e ^ is a monotonic function of x which smoothly interpolates between 
the two limiting cases of homogeneous spheres with the dielectric constants e± and e, respectively, in the 
background with the dielectric constant Sb- Therefore, for any x the effective dielectric constant efu* for a 
lattice of coated spheres is always smaller than the largest of the effective dielectric constants obtained for 
the two limiting cases of a lattice of homogeneous spheres. Using a transfer matrix method to calculate 
the polarization factor a c for a coated sphere with an arbitrary number of coatings, one can show that 
this restriction on the values of the effective dielectric constant e c e ^ for a lattice of coated spheres holds 
also in the general case. Indeed, let us consider a coated sphere made out of N spherical shells with the 
dielectric constant Sj 7 j = 1, . . . , N. Then 

e% f < max{ -/,,} (10) 

where e 3 e ^ is the effective dielectric constants of a lattice of homogeneous spheres with the dielectric 
constant Sj in the background Therefore, one would naively expect that a coating has mere effect of 
interpolating between the limiting cases of homogeneous spheres. 

Nevertheless, considerations based on the value of e c e ^ can be deceptive. In figure || we show a plot 
of the relative L-gap width Al for a close-packed fee lattice of coated spheres as a function of ri/r s , the 
ratio of the interior and whole sphere radii. The solid line corresponds to the case when the refractive 
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index of the core sphere is n\ — 2 (ZnS) and that of the shell is n = 1.45 (silica). The dashed line shows 
the reversed case, i.e., ri\ = 1.45 < n = 2. Obviously, the cases r\jr s = and 1 correspond to the 
limit of homogeneous spheres. Let us denote by A + the larger relative L-gap width of the two limiting 
homogeneous spheres cases. Then in the first case (ni > n) around r\/r s w 0.68 the relative L-gap width 
can be increased by as much as 50% with respect to A + . In the second case (rti < n), the relative L-gap 
width remains smaller than A + for all values of r±/r s . 

In the following, we tested the Maxwell-Garnett formula (Q) against the exact value of the effective 
dielectric constant e c e ^ (refractive index ro e // = ^/e^) obtained directly from a band structure. In Table II 
we collected the effective refractive indices n e ff and n^p (the Maxwell-Garnett value) for those values of 
ri/r s which are nearby the maximum and minimum of the relative L-gap width (cf. figure ^|). Here 
one expects the largest deviations of n^j^ from n e ff because n^jp is a monotonic function of r\jr s and 
hence, it does not have local minima and maxima. 

TABLE II. The effective refractive indices n e // and n^j^ for a simple close-packed fee lattice of coated 
spheres in air. First three columns are for n\ — 2 and n — 1.45, the last three colums are for the 

reversed case n± = 1.45 and n = 2. 



n/r a 


0.65 


0.68 


0.72 


0.75 


0.8 


0.85 


n eff 


1.419 


1.434 


1.444 


1.540 


1.507 


1.469 


MG 
U eff 


1.414 


1.428 


1.438 


1.525 


1.494 


1.460 



It is interesting to note that continues to approximate the exact value n e ff very well. Similarly as 

for the case of homogeneous spheres in air |^8| , in the case of coated spheres in air the Maxwell-Garnett 
formula (Q) slightly underestimates the exact value of the effective refractive index n e ff obtained from a 
band structure 

<f G < n eff . (11) 

4 Conclusions 

We showed that the photonic KKR method is a viable alternative to calculate the photonic band structure 
by reproducing the main features of the photonic band structure obtained by the plane- wave method for a 
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simple three-dimensional fee lattice of homogeneous spheres |8j, |9|, |10[ . We confirmed that for a sufficiently 
high dielectric contrast a full gap opens between the eights and ninth bands, a direct gap opens between 
the fifth and sixth bands, and no gap opens in the spectrum if e& < s s . To obtain a good convergence in 
the frequency range considered, it was sufficient to retain multipole fields with the angular momentum 
up to l m ax — 6. In general, the higher frequency the higher value of l max is needed. In order to reproduce 
the first band and the linear part of the spectrum, l rnax = 1 is enough. The size of a secular equation 
is reduced by almost a factor 10 compared the plane-wave method J6||- Jl0|| which customarily requires 
well above thousand plane waves. Precision of the elements of the secular equation is determined by the 
standard Ewald summation |l(| which yields structure constants p| up to 6 digits. 

For close-packed air spheres, the lower and upper bounds of the full gap are set at the W and X 
points, respectively. For filling fractions less than 0.7 and close to the threshold dielectric contrast for 
which a full gap opens in the spectrum, the gap width is completely determined at the W point. The 
lowest dielectric contrast 8 = £b/e s for which a full gap opens in the spectrum is found to be 8.13 which 
occurs for a close-packed fee lattice. This value is slightly lower than 8.4 obtained by the plane-wave 
method Jl(| . The maximal value of the relative gap width approaches 14% in the close-packed case and 
decreases monotonically as the filling fraction decreases. For readers interested in the the role of Mie's 
resonances in the formation of a band structure we refer to [ j3l| . 

An open question remains what causes such a different behaviour of a lattice of dense spheres compared 
to that of air spheres. This can be partially attributed to the fact that, for a lattice of spheres, £& no 
longer describes the dielectric constant of the surrounding medium, which is instead described by the 
effective dielectric constant £ e ff- Therefore, the bare dielectric contrast 8 is renormalized to 8 e ff = 
max.(s s /s e ff,£ e ft/e s ), where 1 < s e ff < 8 for e s ^ £b- Given the bare dielectric contrast 8, one finds that 
the renormalized dielectric contrast 8^ = e s /e e ff in the case of dense spheres is always smaller than the 
renormalized dielectric contrast 8^ = e e ff in the case of air spheres (see figure ||). The latter is easy 
to verify in the limit of very high bare dielectric contrast 8 — > oo, where Maxwell-Garnett's formula (|7]) 
implies 

8% f ~ 8 (1 - /)/(l + 2/) < 8l n ~ 8 (1 - /)/(l + J/2). 

Nevertheless, a full understanding of the differences between the lattices of air and dense spheres still 
remains a theoretical challenge. 

Such as in the case of homogeneous spheres in air, also in the case of a simple fee lattice of coated 
spheres in air no full gap opens in the spectrum. However, we showed that already a single coating 
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can enhance the relative L-gap width by as much as 50% with respect to A + , the larger relative 
L-gap width of the two limiting cases of homogeneous spheres. This suggests the use of coated spheres 
to reduce the threshold dielectric contrast for which a full gap opens in the spectrum in the case of the 
so-called complex fee lattice (having more than one scatterer in the primitive cell), which will be dealt 
with elsewhere. Note that in contrast to a simple fee lattice, in the latter case a full gap does open in the 
spectrum [To| ]. In general, A^ is not a monotonic function of the effective refractive index n e // and, for 
a single coating, shows one local maximum (minimum) as r\/r s is varied from zero to one (cf. figure ^). 

It is interesting to note that the Maxwell-Garnett formula (Q) provides a very good approximation 
to the exact value of the effective refractive index n e // also for coated spheres. Similarly to the case of 
homogeneous spheres in air pq ], in the case of coated spheres in air the Maxwell-Garnett formula p7[ | 
slightly underestimates the exact value of the effective refractive index n e ff obtained directly from a band 
structure. The Maxwell-Garnett formula holds well beyond the limit Lur s <C 1, r s being the sphere radius, 
for which it was originaly derived. It describes n e ff well almost up to the lowest L gap and can also be used 
to determine the L-midgap angular frequency uj c |3C| . In the case of coated spheres considered here (cf. 
figure ||), the ratio n e ffUJ c /kL £ (0.973,0.995), where fc^ is the length of the Bloch vector corresponding 
to the L point. 
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Figure 1: Photonic band structure for a close-packed fee lattice of air spheres in a background dielectric 
medium with Sb = 9 (n^ = 3). Frequency is plotted in dimcnsionlcss units, A is the lattice constant, c 
is the speed of light in the vacuum. Only a single gap with the central gap frequency v = 2.796 and the 
width Av = 0.044 opens in the spectrum. 



14 



The DOS for a close— packed fee lattice 




0.0 0.5 1.0 1.5 2.0 2.5 
Frequency [c/7tA] 



3.0 3.5 



Figure 2: The DOS per primitive cell for a close-packed fee lattice of air spheres in a background dielectric 
medium with £6 = 9 (rib = 3). Note the gap in the spectrum centered at v = 2.796. 
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Figure 3: Photonic band structure for a close-packed fee lattice of air spheres in a background dielectric 
medium with Eb = 2.1609 (rib = 1.47) - experimental setup reported in [11]. 
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Figure 4: Relative band-gap width, which is the band gap width divided by the midgap frequency, as a 
function of the refractive index contrast for different filling fractions /. 
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Figure 5: The relative L-gap width for a simple closed packed fee lattice of coated spheres as a function 
of the ratio r\jr s of the interior and whole sphere radii. The solid line corresponds to the case when the 
refractive index of the core sphere is ri\ = 2 (ZnS) and that of the shell is n = 1.45 (silica). The dashed 
line shows the reversed case, i.e., n\ = 1.45 and n = 2. Note that in the first case the relative L-gap 
width can be increased by as much as 50% if r\/r s ~ 0.68. The cases r\/r s = and 1 correspond to the 
limiting case of homogeneous spheres. 
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Figure 6: In a medium with finite density of scattering spheres, such as for spheres on a lattice, the bare 
refractive index contrast gets renormalized. This may give a partial explanation of the fact why the band 
structure of an fee lattice of the dense and air spheres shows such a different behaviour. This figure shows 
a typical behaviour of the effective refractive index contrast as a function of the bare refractive index 
contrast for air spheres in a dielectric (solid line) and dense dielectric spheres in air (dashed line). Note 
that the latter is always smaller than the former. This plot was made for filling fraction / = 0.4. The 
Maxwell-Garnett overestimates e e ff for the case of air spheres and underestimates s e ff for the case of 
dielectric spheres in air. As the result, the exact curves are rotated slightly to the right with respect to 
the case when the Maxwell-Garnett value is taken for e e ff. 
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